The Entanglement of Formation for Isotropic States 
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We give an explicit expression for the entanglement of formation for isotropic density matrices 
in arbitrary dimensions in terms of the convex hull of a simple function. For two qutrit isotropic 
states we determine the convex hull and we have strong evidence for its exact form for arbitrary 
dimension. Unlike for two qubits, the entanglement of formation for two qutrits or more is found 
to be a nonanalytic function of the maximally entangled fraction in the regime where the density 
matrix is entangled. 



One of the main goals in quantum information theory 
is to develop a theory of entanglement. A cornerstone of 
this theory will be a good measure of bipartite entangle- 
ment. Such a measure must obey the essential property 
that the entanglement of a bipartite density matrix p 
which is shared by Alice and Bob cannot increase, on 
average, under local quantum operations and classical 
communication (LO + CC) between Alice and Bob. In 
this way, the entanglement captures the truly quantum 
correlations in a bipartite density matrix. For pure bi- 
partite states a good measure of entanglement has been 
found, it is the following quantity: 



E{\^\)= S(Tr B (\^\)), 



(1) 



where S(p) is the von Neumann entropy of p, i.e. S(p) = 
— Tr p logp and Trs(|^) (ip\) is the reduced density ma- 
trix that we obtain by tracing out over Bob's quan- 
tum system. This measure E is unique QJ^] if one re- 
quires the entanglement to obey a set of natural prop- 
erties, such as convexity, non-increase under local mea- 
surements, asymptotic continuity, partial additivity and 
normalization. Moreover, E is a measure of the asymp- 
totic entanglement costs ||| of making the state out of 
a canonical set of states, which we can choose to be EPR 
singlets ^(|01) — 1 10}), which have E = 1. This process 
is reversible, in the sense that one can concentrate || a 
set of n states 1^) with entanglement E to a smaller set 
to = En EPR singlets. 

The situation for mixed states is much more complex. 
In Ref. Q a first measure of mixed state entanglement, 
called the entanglement of formation, was introduced. 
This measure is a candidate for measuring the asymp- 
totic costs of making the density matrix out of a supply 
of EPR singlets. There are no mixed density matrices for 
which this statement has been proved, but neither have 
counterexamples been found so far. The search for a pos- 
sible discrepancy between the entanglement of formation 
and the asymptotic entanglement costs is hampered by 
the fact that we know the entanglement of formation only 
for two qubit systems; Wootters M found an analytic ex- 



pression for the entanglement of formation for all two 
qubit density matrices. 

In this Letter we present the first calculation of the en- 
tanglement of formation of a class of density matrices in 
dimensions higher than C 2 ®C 2 . We explicitly determine 
the entanglement of formation for two qutrit density ma- 
trices in this class and we find an expression in arbitrary 
dimension in terms of the convex hull of a simple func- 
tion. We conjecture the explicit form of this convex hull, 
which can be easily verified in a given dimension. Sur- 
prisingly, the entanglement of formation is found to be a 
nonanalytic function of the parameter characterizing the 
class of states that we consider. 

Let us start by recalling the definition of the entan- 
glement of formation. Let £ p = {pi, \ipi)} be an en- 
semble of pure states which form a decomposition of 
p = 'J2 i Pi\ipi)(4'i\- The entanglement of formation for 
mixed states p is defined as 



E( P ) 



e={pi,\1>i)} ■ 



iPi E(\^)(<Pi\)- 



(2) 



In this Letter we will consider the class of density ma- 
trices, sometimes called isotropic density matrices, which 
are convex mixtures of a maximally entangled state and 
the maximally mixed state: 



PF 



1 - F , 



(3) 



for < F < 1 and |*+) = E<=i R- For F < 1 / d 
these density matrices are separable ||. The entangle- 
ment of formation for states with d = 2 is equal to |I| 



= l + VF(l-F), 



(4) 



where H 2 (.) is the binary entropy function. The states pp 
have the important property || that they are invariant 
under the operation U ®U* for any unitary transforma- 
tion U. The LO + CC "twirling" superoperator S u ® u * 
is defined as 



1 



s 



(P) 



1 



7o/({7) 



dUU ®U*pU ] ®U* 



(5) 



In Ref. 0] the Schmidt number of the isotropic states 
was determined. Instead of making an isotropic state 
out of a set of maximally entangled states, we ask how to 
construct an isotropic state with a given F out of some 
state characterized by a Schmidt vector p. So, let us 
take an arbitrary initial pure state \ip) = Yli=i \[Wi\ a %i &<) 
and consider the effect of twirling. We can write = 
Ua ® Ub J2i yfJn\h i) an d thus 

s V9tr \ Y\1^' ! «■■''■ «,.'>, 

s u ® u * ((i ® v) ® v*)), (6) 

where V = U\\Jb- We define v% = (i\ V \j). The twirled 
state becomes 



-p+ 



i-IE^I 2 /^ 
d 2 - i 



0--P+), 

(7) 



where Wj = ^fplvn and P + = |\f r+ )(\f r+ |. When we choose 
V = lwe find the density matrix pf at F = yfpl] 2 /d. 
For general one can bound 



£^ 2 <[£h] 2 <[£v^] 2 



(8) 



since \vu\ < 1 for all i. Thus the largest value for F is 
obtained by choosing the initial state Ei=i \fjH\H)- 

The use of symmetry makes it possible to give a sim- 
plified expression for the entanglement of formation for 
isotropic states: 

Lemma 1 The entanglement of formation for isotropic 
states in C d <£> C d (d > 2) for F G (1/d, 1] is given by 



E(p F ) = co(i?(F)), 



(9) 



where co{g) denotes the convex hull of the function g and 
R{F) is defined as 



fl(F)=minj#(/2)|F=[V. 



(10) 



and [2 is a Schmidt vector. 

Proof Assume that there exists an optimal decomposition 
of Pf formed by the ensemble {pi,\'4>i(jT'))}, where /P 
denotes the Schmidt vector of the state By twirling 
the l.h.s. and r.h.s. of the equation pp = ^2 Pi\ipi) 
we obtain that pp = J^iPiPFi = z 9 ^ P f where 



T,k V kkVP\ 



> 1 2 



as in Eq. (f7|). Since the decomposition is optimal, each 
Schmidt vector ft 1 has minimal entropy under this con- 
straint. Consider the function 



Rv(F) 



mm{H(p) | F 

P 



d 
i=l 



/p-\ 2 /d}. (12) 



An optimal decomposition of pf is a convex combina- 
tion of pure states each of which corresponds to a certain 
F under twirling. Thus the entanglement of formation 
E(pp) can be obtained by taking the convex hull of the 
functions co{Rv(F)). We can make an additional simpli- 
fication. Eq. (|) implies that R V {F) = Ri(F') = R(F') 
where F' > F for every V. Thus instead of taking the 
convex hull of all functions co(Rv(F)), we can take the 
convex hull of function R + (F) = min x {_R(x) | x > F}. 
In Lemma 2 we will determine R(F) and it is not hard 
to show that R(F) is a monotonically increasing func- 
tion of F. It follows then that R+(F) = R(F) and 
co(R+(F) = co(R(F)). □ 

We now determine the function R(F) defined in Eq. 
(|Io|). Since all the equations are symmetric in pi, we can 
restrict ourselves to solutions which satisfy p± > pi > 
• • • ^ Pd- With the method of Lagrange multipliers we 
get a necessary condition for the minimum 



- 1 - log Pi 



At -I 

Ai + yM, 2 =0, 



(13) 



where Ai, A2 denote the Lagrange multipliers. For fixed 
Ai, A2 this determines the whole set {pi}. Setting pi = 
-ij we obtain an expression of the form log qi — Aqi + 
B where A, B only depend on A!,A 2 . Since a convex 
and a concave function cross each other in at most two 
points, this equation has maximally two possible nonzero 
solutions for qi. Therefore all Schmidt vectors p that 
are possible candidates for the minimum have to satisfy 
the condition pi g {7,(5,0}. Let n be the number of 
entries where pi = 7 and m the number of entries where 
Pi = S. The minimization problem has been reduced 
considerably: For fixed n, m, n + m < d, we minimize 
the function 

nh(i) + mh(S) , (14) 

where h(x) = — x log x, under the constraints 

{717 + mS = 1 , n-y/7 + mVS = V dF } , (15) 

The constraints give rise to a quadratic equation in y/7 
which provides two possible solutions for 7 for every 
choice of n,m: 



(F) 



QiF n ± y mn(m + n — dF) 
n(n + m) 



(16) 



With the first constraint we get the corresponding 



(H) <5± m (F) = (1 - n 7 ± m (F))/m. Since 7 



2 



function in Eq. (jlj) takes the same value for 7+ m and 



Therefore we can restrict ourselves to the solutions 
= Jnm- The pointwise minimum over all possible 



Tmn 

choices for n.mof 



Rnm(F) = H 2 {nq nm ) + ny nm log 1- logm, (17) 

m 



defined on the domain % 

a 



< F < ^^f 1 , is the required 
function R(F). The restriction on the domain comes 
from requiring that 7„ m is a proper solution of Eq. Jl6| ) 
which implies that F < n+ d rn . On the other hand we 
demand that S nm > which implies that F > ^. In this 
regime one can verify that 7„ m > S nm . 

When m = 0, 7 and F are uniquely determined by the 
constraints, i.e. F = Since Rno(j) = Rn'm'(^) for all 
n' + m! = n, we can neglect these cases. 
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FIG. 1. The Rnm(F) functions for d = 3. The solid line 
is Ri2{F), which is the minimal one. The dotted and dashed 
lines are Rn(F) and R2i(F) respectively. 



When d — 3, what remains is a minimization over the 
three functions Ri2(F), R2i{F) and Ru(F), which are 
plotted in Fig. |. For d = 3 we get R(F) = R 12 {F) §. 
Thus the optimal vector fl is always of the form 



(18) 



satisfying 7 > S. 

The case d = 3 is the important one, since it turns out 
that we can relate all the higher dimensional minima to 
d = 3 and prove that 

Lemma 2 For d > 3 the function R{F) = R 1 , d _ 1 (F). 
Proof The case d = 3 is discussed above. Note that 
R\.d-i{^/ d) = 0, which is clearly minimal, so we pro- 
vide a proof for F > 1 / d. Let the minimum be attained 
in d > 3 dimensions by a vector /2 = {/ii}. Let us select 
some subset of the entries of ft, the set {^v }i=i > where 
^2j=i (J-ij — ^ < 1. Since ft is the minimum, it follows 
that the set {Hi^^i is the minimum when we keep the 
other entries of the vector ft fixed. Let [>!■ :— The 
vector yJ is the solution for the minimization of 



h(k), 



(19) 



under the constraints X^=i A*j 
where 



1 and J2 j= i \M = C 



IdF 1 ^ 



(20) 



This last equation can always be written as C = V 'd'F' 
for some F' . Thus the restricted minimization problem 
is equivalent to a d'-dimensional version of the original 
problem, up to the scaling factor k and the additive term 
h(k). When F' < we know that solution of this min- 
imization problem is given by a Schmidt vector fS! which 
corresponds to an unentangled state, i.e. it is of the form 
ft = {1, 0, . . . , 0}. Let us choose three arbitrary ^ out of 
the optimal vector fl. When the resulting F' < |, it fol- 
lows that fl' — {1,0,0}. When F' > |, the three entries 
of fH have to satisfy Eq. (|l8|). So in fact, in both cases 
they satisfy Eq. (|18|). Suppose now that one entry of fl 
is equal to zero. Then it follows that fl cannot have two 
nonzero entries since this would violate condition (^), in 
other words it must be that /I={1,0,...,0}. But this is 
a solution for F = 1 / d. Therefore we get n + m = d for 
F > 1/d. Suppose that n > 2. We can choose the vector 
{7,7,(5} satisfying 7 > 5. Then condition (18) implies 
that 7 = 5. This implies that all entries of ft are identi- 
cal, or/2={i,...,i}. This corresponds to a maximally 
entangled state, which is the unique solution for F = 1. 
Therefore n = 1 and m = d — 1 . □ 

Lemma 1 and Lemma 2 together result in 
Theorem 1 The entanglement of formation E(pf) for 
isotropic states in C d <g> C d (d > 2) for F € {1/d, 1) 
is given by 



E{p F )=co{R 1A -x(F)), 



(21) 



where 



R hd -i(F) = H 2 ( 7 (F)) + (1 - 7(F)) log (d - 1), (22) 
with 

7 (F) = \(VF+ y/(d-l)(l-F)) 2 . (23) 
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FIG. 2. The first and second derivative of Ri,d-i{F) for 
d — 3. The solid line is the second derivative, which is going 
to — oo for F = 1. The dashed line is the first derivative. 



For d — 3 the first and second derivative of the func- 
tion R\2(F) are plotted in Fig. [| The figure shows 
that the function Ri2(F) is not convex near F = 1; its 
second derivative is not positive. In order to determine 
co(i? 12 (F)) for d — 3 we solve the following equations. 
Let Eu ne (F) = aF+log3 — a be the line crossing through 
the point (l,log3). We solve (1) E Une (F) = R 12 (F) and 



(2) = a = f or a an d F. Figure | indicates 

that R\2{F) is monotonically increasing and that there 
is only one region where R\2{F) is not convex, namely 
near F = 1. Therefore the solution to the equations will 
be unique: we find that F = 8/9 and a = 3. For higher 
dimensions, we conjecture, based on examining these two 
equations, that the entanglement of formation in C d ®C d 
is given by 



E(pf) = { 



0. 



#M-t(n 

«^(F-l)+logd, Fe[^l,l] 



(24) 



The correctness of this solution can easily be verified 
for a given d by plotting the function Ri t d-i(F) and its 
second derivative and noting the convex hull of Ri,d-i(F) 
is obtained by calculating where Ri.d~i{F) meets the 
line going through the point (F = 1 , E = log d) and the 
tangent of R[ d -i(F) equals the slope of this line. 

It is surprising to find that E(pp) is nonanalytic in 
the region where pf is an entangled density matrix. An- 
other feature of our solution is that for, say, d = 3 and 
F > 8/9 an optimal decomposition of pf is not one in 
which every pure state has an equal amount of entangle- 
ment. Indeed, the optimal decomposition that we find, is 
a mixture of the maximally entangled state and the en- 
semble of states \ip) obtained by twirling, each of which 
has entanglement E = — 1/3 + log 3. Since every state 
in the optimal decomposition of pf has, under twirling, 
a value of entanglement on Ri_d-i(F), every optimal de- 
composition of pi? for d = 3 in the range F > 8/9 will 
be a mixture of the maximally entangled state and some 
less entangled states. This is in contrast with optimal 
decomposition for E for two qubits. For F > 8/9 more 
than d 2 — 9 pure states must be used in the optimal 
decomposition of pp. We make pf from a maximally en- 
tangled state and the state Pf=8/9 which has rank 9, and 
thus needs at least 9 states in its optimal decomposition. 
In total, this gives 10 states. For F > 8/9 there is no 
optimal decomposition with fewer states: one always has 
to mix in the maximally entangled state with some prob- 
ability. The remaining state p' F either has rank 9 (like 



Pf=8/9) or a lower rank. If it has a lower rank, it must 
be separable, which would imply that the optimal de- 
composition is made from mixing a separable state with 
a maximally entangled state which we know to be false. 
This is the first example of an entangled state for which 
it is proved that the number of pure states in the optimal 
decomposition exceeds the rank of the state (see Ref. |l(]] 
for separable states with this property). 

Crucial in our method is the invariance of the isotropic 
states under a symmetry group of local operations. A 
result similar to Lemma 2 will hold for example for the 
class of Werner states || which are invariant under the 
transformation U <g> U for all U G U(d). Let S u ® u be 
defined as in Eq. (|J), but with omission of the complex 
conjugation. The Werner states pf are characterized 
by a single parameter p. One can prove that E{pf) = 
co(R(p)) where 

R(p) = mm{E(W(i;\) \ S u ® u (\ip)&\) = pf} . (25) 

It may thus be possible to carry out a similar analysis as 
was done here for the Werner states. In a further general- 
ization one could consider the entanglement of formation 
for g ® g or g ® g* invariant states where g G G and G is 
a subgroup of U(d). 
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